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Abstract 

We consider self-similar approximations of non-linear hyperbolic sys- 
tems in one space dimension with Riemann initial data, especially the sys- 
tem dtu^ + A{uir) OxU^r — etdx{B{us) d^u^), with e > 0. We assume that 
the matrix A{u) is strictly hyperbolic and that the diffusion matrix sat- 
isfies \B{u) — ld| << 1. No genuine non-linearity assumption is required. 
We show the existence of a smooth, self-similar solution = u^{x/t) 
which has bounded total variation, uniformly in the diffusion parameter 
e > 0. In the limit e ^ 0, the functions converge towards a solution 
of the Riemann problem associated with the hyperbolic system. A simi- 
lar result is established for the relaxation approximation dtu'^ -{-dxV^ ~ 0, 
dtv'' + a'^ Biu)d:,u'' = t). We also cover the boundary-value 

problem in a half-space for the same regularizations. 

1 Introduction 

In this note we present a continuation of the work by the authors [5]-[5] on self- 
similar regularizations of the Riemann problem associated with a non-hnear 
strictly hyperbolic system in one-space dimension: 

dtu + A{u)d.j;U^O, u = u{t,x) eBs„, t>0,xeIR, (1.1) 

with piecewise constant, initial data 

m(0, x) = ui for x < 0; Ur for a; > 0, (1-2) 

where ui^Ur are constant states in Bs^. Here, Bs^^ C iR^ denotes the ball 
centered at the origin and with radius (5o > 0, and, for all u e Bs^, A{u) is 
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assumed to have distinct and real eigenvalues Ai(u) < . . . < AAr(u) and basis of 
left- and right-eigenvectors lj{u),rj{u), (1 < j < N). 

Following Dafermos [1] [5] and Slemrod [T3] who advocate the use of self- 
similar regularizations to capture the whole wave fan structure of the Riemann 
problem, we consider solutions constructed by self-similar vanishing diffusion 
associated with a general diffusion matrix B = B{u), that is we search for 
solutions of 

dtu^ + A{u^)dxU^ — etdx{B{u^)dxU^'), e > 0. (1.3) 

Due to the choice of the scaling et, this system admits solutions — u^{x/t)^ 
and, therefore, we refer to (|1.2p - (|1.3p as the self-similar diffusive Riemann prob- 
lem. The matrix B = B{u) is assumed to depend smoothly upon u and to 
remain sufficient close to the identity matrix, that is, for some given matrix 
norm and for rj > sufficiently small 

sup \B{u) - Id| < r]. (1.4) 

The method of analysis introduced below is not a priori restricted to (|1.3p - (|1.4[) . 
and generalizations are discussed at the end of this note and in [51 H] . 

The techniques developed so far for general hyperbolic systems (see [TH [T^ 
and [5]-[7j) were restricted to regularizations based on the identity diffusion 
matrix. The new approach introduced here allows us to cover classes of approx- 
imations based on general diffusion matrices (or relaxation terms, see below). 
This degree of generality is especially important for non-conservative systems 
[TUl fTT] and for the boundary- value problem [5] , whose solutions are known to 
strongly depend upon the specific regularization. 

2 Main results 

By the property of propagation at finite speed, a self-similar solution u = u(y) 
of the Riemann problem is constant outside a sufficiently large, compact interval 
[— L,L], i.e.: u{y) = ui for y < —L and u{y) — Ur for y > L. As is customary, 
we assume that So is sufficiently small so that the wave speeds Xj{u) remain 
close to the constant speeds Xj{0) and are uniformly separated in the sense that 

Aj<Aj(m)<Aj, ueBso, 

for some constants —L < Aj^ < Ai < A2 < . . . < A^ < Ajv < L. 

We establish the following theorem. 

Theorem 2.1. Consider the non-linear, strictly hyperbolic system il.l]) to- 
gether with its parabolic regularization \1.3\) - [L4^ . There exist (sufficiently 
small) constants Si,ri > and a (sufficiently large) constant Co > such 
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that for any initial data ui,Ur £ Bg^ the self-similar diffusive Riemann prob- 
lem lil.2\) - nT^) admits a smooth solution — u^(x/t) S Bs^ defined for all 
y — x/t G [— L, L], which has uniformly hounded variation, 

TV\{u,)<Ca\ur-ui\, 

and converges strongly to some limit u : [~L,L] — > Bsg-' 

u!^ ^ u in the norm, as e 0. 

The limit function satisfies the following properties. The function y t— > u{y) has 
bounded total variation, that is, TV:'^]^{u) < Cq \ur — ui\, and is constant on each 
interval [Aj,Aj^i]. If U.l\) is a system of conservation laws, i.e. A = Df for 
some flux f : Bsg , then the limit is a distributional solution of 

dtu + d^fiu) ^0. (2.1) 

// U, F) : Bs^ IR X IR^ is an entropy / entropy flux pair associated with V2.1]) 
and the diffusion matrix satisfies the convexity-like condition V^U • B > 0, then 
the solution u satisfies the entropy inequality 

dtU{u) +d^F{u) <0. (2.2) 
We have also the following description of the wave curves. 

Theorem 2.2. With the notation and assumptions in Theorem \2.1\ to each 
j- characteristic family and each left-hand state ui one can associate a j-wave 
curve 

yVj{ui) {ur = ■>jjj{m; ui) / m e {mj,rnj)}, 

issuing from ui, which, by definition, is made of all right-hand states Ur at- 
tainable by a Riemann solution u = u{y), with left-hand state ui, by using only 
j -waves, that is such that 

u{y) ^ ui for y < kj] Urfory>kj. 

Moreover, the mapping ipj : (rrij , ruj ) x Bs^ — *■ Bs^ is Lipschitz continuous with 
respect to both arguments, and for some small constant c > 

dm'4>3{rn;ui) £ Cj := {w e IR^ / \w lj{{))\ > (1 - c) 

Moreover, the characteristic component y ^ ctjiu) ■— Iji^) ' u'{y) is a non- 
negative measure in the interval [A^ , A^]. 

For the proof of these results as well as a characterization of the limit when 
(jl.ip is a general non-conservative system, we refer to [8j |9] . 
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We will here sketch the proof of Theorem 12.11 To handle general diffusion 
matrix B{u), the following generalized eigenvalue problem is introduced 

{-y + M^)) rj{u, y) = y) B{u) rj{u, y), 

hiu,y) ■ {~y + A{u)) = fJ.j{u,y) lj{u,y) ■ B{u). 

In view of (|1.4p . one has 9j (u. y] — rj{u) + 0{ri) and lj{u,y) = lj{u) + 0{ri). The 
proof relies on a suitable asymptotic expansion of the solution = Ug{x/t), of 
the form 

u'^ — a'^jrj{ue, •) with := lj{ue, •) • u'^. 

3 

Omitting e, we deduce that the components satisfy a coupled system of N 
differential equations: 

a'i - ^ + X! ■) = Qiiu; •) := ^ Kijk{u, ■) Qj ak, 

where 

TTy (w, ■) := li{u, •) • B{u) dyrj{u, ■), 
Kijk{u,-) := ~li{u,-) ■ Du{Brk){u,-)rj{u,-). 
The system under study has the form 

4 - a. + OM J2 = 0(1) J2 

In a central part of our argument we study the homogeneous system 







^TTij{u,-)(pj ^0, ifi^ {ipi,...,ipN), (2.3) 



and establish that it has solutions ipj, referred to as the linearized i-wave mea- 
sures associated with the function u, which are "close" (in a sense to be specified) 
to the following normalized solutions of the corresponding uncoupled system 
(obtain by setting ry = 0) 

g-ffi/e /"^ _ / r 

ip* -.^ — J — , h -^ e^'/'^dy, gdy) ■■= - fi^{u{x),x)dx. 

J^i J-L J Pi 

Here, the constants pi are determined so that the functions gi are non-negative. 



Theorem 2.3. The system 112. 3\) admits a solution if such that for all i = 
1, . . . , TV and y G [—L, L] 

(l-0(r;)) ^Uy)-eO{rj) J^^^iy) < ^.{y) < (l+0(r;)) ipUy)+eO{ri) ^^*{y). 
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In constrast with the functions the functions ipi need not be positive. 
Next, to control the total variation of the solutions of (|1.3p . we derive Glimm-like 
estimates on the wave interaction coefEcients 



for some constants Ci € [Aj,Ai]. We gain useful information on the possible 
growth of the total variation of solutions. Roughly speaking, the coefficient 
Fijk bounds the contribution to the i-th family due to interactions between 
waves of the j-th and k-th characteristic families. 

3 Generalizations 

The results above have been also extended to relaxation approximations and 
boundary-value problems. In particular, we can handle relaxation approxima- 
tions associated with the conservative system p.ip 



where — u'^{x, t) and — v'^{x, t) are the unknowns, and e > is a relaxation 
parameter. 

We also study (|1.3p and (|3.ip in the presence of a boundary, when there exists 
an index p such that < Ap, and that at most one wave family is characteristic, 
that is, G (Ap,Ap). We consider (jl.ip on the interval y G [0, L], and prove 
the existence of a solution with uniformly bounded variation. To handle the 
boundary layer, we modify the previous definition of the functions V'^, j < P 
and carefully estimate the coefficients F*ji,{y) when e ^ 0. 

In addition, following pioneering work by Fan and Slemrod j4j who studied 
the effect of artificial viscosity terms, we consider a system arising in liquid- 
vapor phase dynamics with physical viscosity and capillarity effects taken into 
account. We establish uniform total variation bounds, allowing us to deduce new 
existence results. Our analysis cover both the hyperbolic and the hyperbolic- 
elliptic regimes and apply to arbitrarily large Riemann data. The proofs rely 
on a new technique of reduction to two coupled scalar equations associated 
with the two wave fans of the system. Strong convergence to a weak solu- 
tion of bounded variation is established in the hyperbolic regime, while in the 
hyperbolic-elliptic regime a stationary singularity near the axis separating the 
two wave fans, or more generally an almost-stationary oscillating wave pattern 
(of thickness depending upon the capillarity-viscosity ratio) are observed which 
prevent the solution to have globally bounded variation. 




dtu' + d,v' = 0, 



ifiu^-v^), (3.1) 
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